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Abstract: Solar simulators are common laboratory devices to artificially reproduce the Sun emission spectrum. Their use in 
optic tests allows to study the effect of the solar radiation on both materials and components. This paper focuses on the 
effective design of the ellipsoidal reflector for concentrating solar simulators. A ray -tracing analytic model integrated to 
Monte Carlo simulation is proposed as an effective approach to optimize the reflector geometric configuration and to 
maximize the target incident radiation level and flux distribution. The developed model reproduces the ray trajectories from 
the source to the target and it includes the physical and optic phenomena affecting the light rays, e.g. absorption, deviation, 
reflection, distortion, etc. A realistic case study for the effective design of the ellipsoidal reflector for a single source small 
scale simulator, integrating a commercial xenon short arc lamp, is proposed to both validate and apply the proposed 
approach. Several scenarios are tested and the main obtained evidences are summarized. 
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I. Introduction 

Solar simulators provide a luminous flux approximating natural sunlight spectrum. Their basic structure includes a 
metal support frame, a light source, e.g. high flux arc lamp with power supplier, igniter and electronic load, and a reflective 
surface to properly orient the emitted rays lighting the target area. The reflector shape allows the system to generate a 
concentrated or non -concentrated light beam through an ellipsoidal or parabolic reflector. Fig. 1 shows a 3D scheme of a 
concentrating solar simulator, highlighting its functional modules. 
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Figure 1: Single Source Solar Simulator, Functional Module*. 



The relevance of such systems, in laboratory tests and analyses, is frequently discussed by the recent literature 
presenting several applications for a wide set of research fields. Petrash et al. [1] both review the topic and describe a 1 1000 
suns high-flux solar simulator. Domfnguez et al. [2, 3], Pravettoni et al. [4], Rehn and Hartwig [5], Hussain et al. [6] and 
Meng et al. [7] present different studies about the design and development of high flux solar simulators applied to both 
concentrating and non -concentrating photovoltaic systems. 

Amoh [8] and Meng et al. [9] describe the design of solar simulators to test multi-junction solar cells for terrestrial 
and space applications. Kreuger et al. [10] develop a 45kW solar simulator for high-temperature solar thermal and thermo- 
chemical researches, while Codd et al. [11] present a low cost high flux simulator to study the optical melting and light 
absorption behavior of molten salts. All contributions focus on the relevance of the proper design of the system to achieve 
high performances in both flux intensity and uniformity on the target area. The mirror reflective surface represents a crucial 
component to gain these purposes. An accurate shape design and simulation of the physical and optic properties is essential 
before the simulator construction [12]. 

Ray-tracing algorithms are recognized as effective approaches to test the performances of different configurations 
of the reflective surface [1, 13-15]. 

This paper presents a ray-tracing model integrated to a Monte Carlo simulation to effectively design the ellipsoidal 
reflector of a single source concentrating solar simulator. The proposed model analytically reproduces the ray trajectories 
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from the source to the target and it includes the physical and optic phenomena and the distortions affecting the light rays. 
Furthermore, the ray-tracing model integration to the Monte Carlo simulation allows to test and compare several 
configurations of the reflector geometry and to study the most performing shape. This paper discusses such an approach and 
applies it to a realistic case study. The design of the reflector for a small scale solar simulator based on a commercial xenon 
short arc lamp is assessed. Particularly, the description of the implemented steps is provided together with the analysis of the 
simulated scenarios. Finally, the major outcomes are presented and properly commented. 

The proposed approach goal is to support the designer of solar simulator optic systems with an integrated and easy- 
use method to rapidly test a system geometry before the prototyping phase and laboratory validation. An integrated approach 
to gain this purpose is of strong interest for both scientists and practitioners and its development is frequently encouraged. 

The reminder of the present paper is organized as follows: Section II describes the steps of the developed ray- 
tracing model, while Section III introduces the Monte Carlo analysis and the aforementioned case study to design the 
ellipsoidal reflector of a single source concentrating solar simulator. The obtained results are presented and widely discussed 
in Section IV. Finally, Section V concludes this paper providing suggestions for further research. 

II. Ray-Tracing Model 

In geometrical optics, the foci of an ellipsoid of revolution are conjugate points [1]. If no distortion effects occur, 
each ray emitted by a punctiform source located in one of the foci passes through the other after a single specular reflection. 
According to this principle, the concentrating solar simulators are designed. The light source, i.e. the yellow spot in Fig. 2, 
reproducing the Sun emission spectrum and the target surface, i.e. the blue spot in Fig. 2, are located in correspondence of 
the two foci of an ellipsoidal surface, while the reflector lies on a portion of such an area and it is limited by a truncation 
section as depicted in Fig. 2 presenting the overall geometry of the system. 



v 



'7 



Figure 2: Single Source Solar Simulator Geometry and Notations. 

Considering the operative contexts, the following conditions and phenomena contribute to reduce the global system 
radiation transfer efficiency, expressed as the ratio between the light flux that reaches the target and the global emitted flux. 
The finite area of the light source; 

The absorption phenomena caused by the presence of the quartz bulb, the electrodes and the reflective surface; 
The deviation and distortion phenomena caused by the specular dispersion errors of the reflective surface; 
The losses caused by the rays falling out of the reflector shape. 

Such conditions cannot be neglected in the solar simulator design. Their impact in reducing the system 
performances is strongly correlated to the features of the emitting source, the target shape and, particularly, to the reflector 
characteristics. From this perspective, a model approaching this issue from an analytic point of view is of crucial interest. 

The proposed ray-tracing approach faces this topic and it studies the ray trajectories, predicting the system global 
performances for a given configuration of the source, the reflector and the target surface. The next Fig. 3 shows the ray- 
tracing model flow-chart. It summarizes the step sequence of the proposed approach together with the stages in which the 
losses in the transfer efficiency generally occur. 
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Selection of the incident ray 



Prediction of the incident ray direction v 
according to Lambert's cosine law distribution 



Calculation of the point of intersection between v 
and the ellipsoidal reflector, P, 

I 

Definition of the reflected ray direction r, 
including reflector distortion effects 

I 



Calculation of the point of intersection between r 
and the target, P 2 



Transmission efficiency and 
flux distribution analysis 



/ Losses at light source stage 
- quartz bulb absorptivity 

- electrode absorptivity and ray deviation 

/ Losses at reflector stage 

- surface absorptivity 

- ray out of reflector shape, no intersec tion 



Losses at target stage 
t of target shape, 



Figure 3: Ray-Tracing Model Flow-Chart. 



According to the major literature [1,2, 10] the light source is assumed to emit isotropic radiation uniformly from its 
surface. Consequently, the emission point, P 0 , is randomly located on the whole source surface. The incident ray direction, v, 
is defined following the Lambert's cosine law distribution, as in (1) [16]. 



>s(si 



(1) 



where n is the normal direction to the emitting surface, in P 0 , and U a random number drawn from a [0, 1] uniform 
distribution. The quartz bulb and electrodes absorption phenomena are considered introducing two coefficients, i.e. the bulb 
and the electrodes absorption coefficients, reducing the emitted rays and decreasing the system efficiency, i.e. the losses at 
the light source stage. 

For each emitted ray, the point of intersection with the ellipsoidal surface, P,, is calculated. If Pi falls out of the 
reflector shape or it falls in the hole, close to the ellipsoid vertex, necessary to install the light source, the ray is considered 
lost and the process finishes. Otherwise, two possibilities occur. Generally, the mirror reflects the ray but, in few cases, an 
absorption phenomenon occurs and the ray is not reflected at all. In such a circumstance, modeled considering a further 
absorption coefficient, the process ends, i.e. the losses at the reflector stage. 

Finally, considering the reflected rays, their direction, r, is calculated. The vector r is affected by the distortion 
effects caused by the specular dispersion errors of the reflective surface. As widely discussed by Cooper and Steinfeld [15], 
the geometric surface errors modify the normal vector, k, to the ellipsoid surface. The authors identify two angular 
components of the dispersion error, i.e. the azimuthal angular component, 8 err , and the circumferential component, (p err . By 
applying the, so called, Rayleigh method they outline the following expressions to estimate them: 



<Pcn 



= 2nU 



(2) 



(3) 



where a err is the standard deviation of the dispersion azimuthal angular error distribution, including all distortion effects, 
and U a random number drawn from a [0,1] uniform distribution. As a consequence, to estimate the direction of k, in the 
point of intersection P u the theoretic normal vector k' needs to be twofold rotated. The rotation angles are 8 err and, then, 
<p err . The former rotation is around a vector orthogonal to the plane where the major ellipse lies, while the latter rotation is 
around k'. 

The prediction of the normal vector to the reflective surface, in P h allows to calculate the reflected ray direction, r, 
according to (4) [16]. 



= v - 2 ■ (k x v) x k 



(4) 



The point of intersection between r and the plane where the target lies allows to find the coordinates of the common 
point P 2 . If P 2 is inside the target area the ray correctly hits the target, otherwise the ray is lost and the transfer efficiency 
decreases, i.e. the losses at the target stage. This study does not consider multiple reflection phenomena. 
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Finally, the distance and mutual position between P 2 and the ellipsoid focus point allows to study the radiative 
incident flux distribution on the target area. 

III. Monte Carlo Simulation 

Several geometric and optic parameters affect the global transfer efficiency of solar simulator systems. A list of them, 
classified according to the physical component they belong to, is provided in the following. 

Light source (generally an high flux arc lamp) 
Shape and dimensions; 
Emission light spectrum; 

Emission surface shape and dimensions, e.g. sphere, cylinder, etc.; 
Emission direction distribution; 

Absorption coefficients of quartz bulb and electrodes (if present); 
Interference angle of electrodes (if present). 

Ellipsoidal reflective surface 

Reflector shape, identified by the two ellipsoid semi-axes or by the major semi-axis and the truncation diameter; 
Reflector length; 
Absorption coefficient; 

Standard deviation of the dispersion azimuthal angular error distribution. 

Target surface 

Shape, e.g. circular, squared, rectangular; 
Dimensions; 

Relative position toward the ellipsoid. 

For a given a set of the previous parameters, the geometric and optical features of the solar simulator are univocally 
identified and the ray-tracing model, described in previous Section II, is applicable to study the global performances of the 
system. In particular, the simulation of a large number of the emitted light rays allows to assess the expected distribution of 
the light on the target and to estimate the system global transfer efficiency. Such a simulation is feasible through the well- 
known Monte Carlo simulative approach when integrated to the described ray-tracing model. 

Finally, varying one or several of the previous parameters, through a multi-scenario analysis, the best configuration 
of the system can be identified. Such a strategy is adopted in the realistic case study described in the following paragraph. 

3.1 Case Study. Design of a Single Source Concentrating Solar Simulator 

The following case study provides an empirical application of the proposed approach. The design of a small-scale 
single source concentrating solar simulator is assessed. The system overall structure is similar to that proposed in previous 
Fig. 1. In particular, the emitting source is a commercial OSRAM XBO® 3000W/HTC OFR xenon short arc lamp with a 
luminous flux of 1300001umen and an average luminance of 85000cd/cm [17]. The target area is assumed squared, centered 
on the ellipsoid focus point and its dimensions are of 50 x 50mm. Such features are assumed constant, while several 
configurations of the reflector shape, corresponding to different sets of the aforementioned parameters, are tested and the 
performances compared. 

Further details about the considered high flux arc lamp are in Table 1 and represented in Fig. 4. The quartz bulb 
absorption coefficient is of 4% and the electrode absorption coefficient is of 98%. 

Table 1: Geometric Features of the Emitting Source. Notations Refer to Fig. 4. 



Lamp length (overall) 


h 


398mm 


Lamp length 


h 


350mm 


Lamp cathode length 


a 


165mm 


Electrode gap (cold) 


CO 


6mm 


Bulb diameter 


d 


60mm 


Electrode interference 




30° 


angles 


h 


20° 
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High Flux Arc Lamp Geometry. 



Considering the ellipsoidal mirror reflector, the next Table 2 summarizes the tested scenarios corresponding to 
different geometric shapes and optic performances. The ranges of variation and the incremental steps for four of the most 
relevant parameters are presented. 

Table 2: Tested Configurations for the Ellipsoidal Reflector [mm]. Refer to Fig. 2 for Notations. 





Min value 


Max value 


Step 


A 


200 


1000 


100 


TD 


100 


2A 


50 


°err 


0.005 


0.01 


0.005 


L 


A-jA 2 -TD 2 /4 


A 


50 



Furthermore, a constant mirror absorption coefficient of 4% is considered in the analysis. 3840 scenarios appear and 
are simulated, i.e. multi-scenario analysis. For each scenario, JV = 5x10 emitted rays are traced and the results collected. 

IV. Results and Discussion 

The measurement of the performances for each simulated scenario is assessed collecting the following data: 
N A , number of the rays absorbed by the light source; 

N H , number of the rays lost due to the presence of the hole used to install the light source; 

N L , number of the rays falling out of the reflector shape; 

N R , number of the rays absorbed by the mirror reflector; 

N T , number of the reflected rays hitting the target; 

JV 0 , number of the reflected rays that do not hit the target. 

Such a data allow to calculate the following key performance indicators highlighting the impact of the reflector 
features on the solar simulator efficiency: 

Losses due to the reflector shape, i.e. the ellipsoid shape and eccentricity, the hole and the truncation diameters; 



Losses due to the optic and distortion effects caused by the reflector surface error distribution; 

f _ N R +N 0 

^2- n _ Na _ Nh _ Nl (6) 

Global reflector transfer efficiency; 

r 7 = (l-f 1 )'(l-f 2 ) = ^ (7) 

Statistical distribution of the reflected rays on the target surface, i.e. the mean distance M D and its standard deviation a B 
between the point of intersection P 2 and the target center. 

The Table 3 shows a subset of the obtained results for the tested scenarios presenting the twenty best and worst 
cases. In addition to the previous notations, e indicates the ellipsoidal reflector eccentricity, as in (8), and included in the 
[0,1] range. 



£ = Jl-B 2 /A 2 (8) 

Furthermore, Fig. 5 shows an example of radiative flux map. It considers the best of the simulated scenarios. The 
squared dashed line identifies the target surface whereas all dots inside the square are the rays hitting the target. The other 
dots are the rays causing the losses at the target stage. 
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Figure 5: Radiative Flux Map for the Best Scenario. 

The values of the global transfer efficiency vary from 92.407% of the best scenario to 1.072% of the worst case. 
Consequently, a first relevant outcome of the analysis is the heavily influence of the reflector design on the solar simulator 
performances. 

Considering the best scenarios, the ^ and f 2 l° ss indices present values lower than 9%, and the large amount of the 
reflected rays are concentrated close to the target, i.e. the mean distance between the rays and the ellipsoid focus point, M D , 
is close to 10mm and the standard deviation, a D , is included between [5,13]mm. On the contrary, the most relevant cause for 
the performance decrease is the reflector shape losses. With reference to the worst scenarios, high values of f ls i.e. greater 
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than 93%, are always experienced, while f 2 does not present a regular trend. Such losses depend on the length of the 
reflector, i.e. the aforementioned parameter L. All the worst scenarios present little values for such a parameter, e.g. 
50-^100mm, so that the majority of the emitted rays is lost because no intersection with the mirror surface occurs. The high 
number of the rays falling out of the reflector shape, N L , included between 80% and 86%, clearly highlights this evidence. 

Furthermore, the standard deviation of the dispersion azimuthal angular error distribution, a err , represents another 
relevant parameter affecting the global performances of the system. As expected, the lower o err , the higher the global 
transfer efficiency values are. The standard deviation error depends on the accuracy of the reflector manufacturing. The 
decrease of this parameter is generally associated to the increase of the reflector production costs. Fig. 6 correlates the 
reflector length to the global transfer efficiency for the two simulated values of o err , i.e. 0.005mm and 0.01mm. 



II II II II II II, 



Reflector length, L [mm] 



Figure 6: Correlation between the Reflector Length and the Global Transfer Efficiency for the Two Values of a t 



The results present similarities in the waveforms. Low values of the reflector length are associated to poor 
performances, i.e. r\ < 30%. Optimal conditions are, respectively, for a reflector length of 800mm and a err = 0.005mm and 
of 600mm for a err = 0.01mm. A significant performance increase occurs for values of L included in [50,450]mm range, 
while for the higher values of the reflector length, i.e. L > 500mm, the global transfer efficiency presents comparable 
values. Finally, considering the gap between the performances in the trends identified by the two values of a err , an increase, 
from 4.241% to 24.466%, occurs. High values of a err have a crucial impact on the global transfer efficiency in presence of 
high values of L. In fact, long reflectors force the emitted rays to sweep out long trajectories from the source to the mirror 
and, then, from the mirror to the target. An anomaly in the trajectories generates an angular deviation of the ray path. Such a 
deviation is amplified by the distance between the mirror and the target. Consequently, if L increases the standard deviation 
of the dispersion azimuthal angular error distribution must have low values not to significantly reduce r\. 

Another relevant parameter affecting the mirror reflector design is the ellipsoid eccentricity, e, defined in previous 
(8) and included in the [0,1] range. It identifies the mutual position of the vertices and the foci. If £ is equal to 0 the ellipsoid 
is a sphere, i.e. A = B. Values of £ between 0 and 1 are for eccentric geometries in which B < A. If £ = 1 the ellipsoid 
degenerates into a plane and the foci lay upon the vertices on the major axis. The developed multi-scenario analysis 
highlights a range of the optimal values for the ellipsoid eccentricity included between 0.75 and 0.9, as in Fig. 7, correlating 
the ellipsoidal mirror eccentricity to the values of r\. Each dot represents one of the 3840 simulated scenarios. 
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Figure Correlation between the Ellipsoid bcccniricil) and Global Transfer Efficiency. 
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This outcome may be in contrast to the major literature evidences suggesting low values of £ to maximize the reflector 
global transfer efficiency [16,18]. On the contrary, in the proposed analysis values of eccentricity close to zero generate the 
worst performances. A reasonable explanation for this evidence lies in the adopted reflector modeling approach. The 
literature ray-tracing models and related results approximate the reflector with an ellipsoid of revolution neglecting both the 
truncation section, i.e. the previously called TD parameter, and the hole necessary to install the light source. The proposed 
ray-tracing approach includes these two elements in the analytical model to provide a realistic and accurate description of the 
physical system. The presence of these elements significantly modifies the geometric and optic features of the solar 
simulator introducing the so-called losses at the reflector stage (see Fig. 3) that significantly contribute to the global transfer 
efficiency decrease, especially for the scenarios in which L and e assume low values (see Table 3). In fact, if £ is low the foci 
are located far from the vertices and close to the ellipsoid geometric center, i.e. the point of intersection of the two axes. In 
this circumstance, the light source juts out from the reflector profile and a large number of the emitted rays does not hit the 
reflector surface. The lower the reflector length, the higher such losses are. 

On the contrary, in eccentric reflectors the light source is close to the ellipsoid major axis vertex and a lower number of 
rays is lost. However, very high values of £, i.e. £ > 0.9, cause a further increase of the losses at the reflector stage and a 
decrease of r]. This is due to the presence of the hole for the light source installation. A focus point located close to the 
reflector vertex increases the value of N H , i.e. the number of the rays lost due to the presence of the hole used to install the 
light source, so that, also in this case, the global transfer efficiency decreases. As introduced before, optimal values for the 
reflector eccentricity are in the [0.75,0.9] range. 

V. Conclusions and Further Research 

This paper presents an approach facing the effective design of solar simulators. The developed model reproduces 
the trajectories of the light rays considering the main physical and optic phenomena that occur from the source to the target 
area. The ellipsoidal reflector geometries are focused. In particular, the solar simulator reflector is a truncated ellipsoid of 
revolution with the light source located on one focus point and the target area on the other. The proposed ray-tracing model 
is integrated to a Monte Carlo simulation to study and compare the performances of several reflector geometries. A case 
study, based on a commercial xenon short arc lamp, is described simulating 3840 scenarios and varying four major 
parameters affecting the quality and reflectivity of the mirror surface, i.e. the ellipsoid major axis, the truncation diameter, 
the reflector length and the standard deviation of the azimuthal angular error distribution. For each simulated scenario, the 
data about the losses and the number of rays on target are collected and summarized in the three key performance indices 
proposed in (6) to (9) together with a statistic analysis of the distribution of rays on the target. 

The main outcomes highlight the relevance of the proper design of the reflector shape to obtain high values of the 
global transfer efficiency. The gap between the best and worst scenarios is higher than 90%. Furthermore, correlations 
between the four considered parameters are highlighted. As example, high values of the reflector length, in presence of high 
values of the standard deviation of the azimuthal angular error distribution, amplify the global transfer efficiency decrease, 
while, low values of the ellipsoid eccentricity cause an increase of the lost rays. 

The obtained parameter values, for the best scenario, are of 600mm for both the ellipsoid major semi-axis and 
reflector length, 325mm for the minor semi-axis and a err equals to 0.005mm. For this scenario, the global transfer efficiency 
is 92.407% while the distribution of rays on the target area presents a mean distance from the ellipsoid focus point of 
8.284mm and a standard deviation of 6.387mm. 

Both the model description and its application to the realistic case study convey to point out the interest of the 
proposed approach for scientists and practitioners. Its application before the system development and manufacture allows to 
simulate the system behavior and to assess the expected performances. Furthermore, the inclusion of the operative geometric 
constraints, e.g. the truncation diameter, the lamp installation hole, etc., to the ray-tracing model increases the approach 
applicability overcoming the assumptions of several of the models proposed by the recent literature. 

Further research deals with a validation of the case study results through the development of the solar simulator and 
a further field-campaign. To this purpose, the authors already purchased the ellipsoidal reflector and they are now 
developing the overall structure of the solar simulator to collect the experimental data to be compared to the evidences 
coming from the proposed approach. At last the enlargement of the approach to multi source solar simulators and the 
inclusion of multiple reflection phenomena is encouraged. 
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